Repeatable light paths in the shearfree normal cosmological models 
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Conditions for the existence of repeatable light paths (RLPs) in the shearfree normal cosmological 
models are investigated. It is found that in the conformally nonflat models the only RLPs are 
radial null geodesies (in the spherical case) and their analogues in the plane- and hyperbolically 
symmetric cases. In the conformally flat Stephani models, there exist special spherically-, plane- 
and hyperbolically symmetric subcases, in which all null geodesies are RLPs. They are slightly 
more general than the Friedmann - Lemaitre - Robertson - Walker models of the corresponding 
symmetries: their curvature index function k{t) and the scale factor R{t) are expressed through a 
single function of time. In addition to that, there exist special cases of the Stephani solution in 
which some of the null geodesies are RLPs. All these special cases are identified. 



I. THE MOTIVATION 

In a recent paper [l[ it was found that in a general 
Szekeres model 0, i] of the 13' ^0 class there are no re- 
peatable light paths (RLPs). This means, given a fixed 
light source S and a fixed observer O, two light rays emit- 
ted from S in such a direction that they hit O, but at dif- 
ferent time instants, intersect different sequences of mat- 
ter world lines on the way. The observer sees then the 
light source slowly drifting through the sky. Depending 
on the mass distribution along the ray, and on the posi- 
tions of O and S with respect to each other, the average 
rate of change of the direction toward S measured by O 
would be between 10~® and 10~^ arc sec per year. With 
the current precision of direction determination equal to 
10~^ arc sec [ij, in the most favourable configuration 10 
years of monitoring would be required to measure this 
effect. Since this drift is strictly zero in the Friedmann - 
Lemaitre - Robertson - Walker (FLRW) models, it can 
be a qualitative observational test of homogeneity of the 
Universe. 

In Ref. [l| it was found that the only subcase of 
the /3' ^ Szekeres family of models in which all null 
geodesies are RLPs are the Friedmann models. The con- 
dition for all null geodesies to be RLPs was the vanishing 
of shear in the flow of the cosmic medium, which reduces 
the Szekeres model to the Friedmann limit. In the gen- 
eral Szekeres models, the matter source moves with zero 
acceleration (i.e. along timelike geodesies) and zero ro- 
tation. This gives rise to the question whether the non- 
RLP phenomenon is caused by shear or by the model 
being non-Friedmannian. 

This question is addressed in the present paper. The 
condition of existence of RLPs is applied here to the 
shearfree normal (SFN) cosmological models [1, [1] - an- 
other class of generalisations of the FLRW models, in 
which the matter source is a perfect fluid moving with 
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zero rotation and zero shear, but nonzero acceleration. 
In the most general conformally nonflat SFN models, 
the only RLPs are radial null geodesies (in the spher- 
ical case) and their analogues in the plane- and hy- 
perbolically symmetric cases. In the conformally flat 
Stephani models, there exist special spherically-, plane-, 
hyperbolically- and axially symmetric subcases, in which 
all null geodesies are RLPs. They are slightly more gen- 
eral than the FLRW models of the corresponding symme- 
tries: their curvature index function k{t) and the scale 
factor R{t) are expressed through a single function of 
time. There also exist special cases of the Stephani 
model, in which some of the null geodesies are RLPs. 
For example, in the general axially symmetric case the 
RLPs are those null geodesies that intersect each space 
t — constant on the axis of symmetry (in full analogy 
with the Szekeres models |lj). All these special cases are 
identified. 

Thus, the key to the RLP property is not vanishing 
shear, but high symmetry. 



II. THE SHEARFREE NORMAL (SFN) 
COSMOLOGICAL MODELS 



The SFN models are solutions of Einstein's equations 
with a perfect fluid source that moves with zero rotation 
and zero shear, but nonzero expansion and acceleration. 
The full collection of these models was first found by 
Barnes [3] in 1973, but special cases were known before 
(see Ref. Q for a classification of special cases and the 
account of historical order) . This family of models con- 
sists of the conformally flat Stephani Universe ^] found 
in 1967, and 3 subfamilies of Petrov type D solutions, 
first found in full generality by Barnes [j] . Each of these 
type D subfamilies has a 3-dimensional symmetry group 
acting on 2-dimensional orbits; the symmetry is either 
spherical or plane or hyperbolic. The spherically sym- 
metric subfamily was first presented by Kustaanheimo 
and Qvist in 1948 0, but one special case of it was de- 
rived by McVittie already in 1933 The other two 
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Petrov type D subfamilies first emerged in tfie paper by 
Barnes Q. This author Q found a coordinate system 
that covers ah 3 subfamilies, but for the present paper it 
will be more convenient to consider them separately. 

In the Petrov type D case, the metric in comoving co- 
ordinates is 



ds^ = 



FV,t 
V 



1 



T72 (da 



dy^ + dzJ 



(2.1) 



where F{t) is an arbitrary function, related to the expan- 
sion scalar 6 hy 9 = 3/F. The Einstein equations reduce 
to the single equation: 



W,uu /W^ = f{u), 



(2.2) 



where f{u) is another arbitrary function, while the vari- 
able u and the function w are related to the coordinates 
X, ?/, 2, and to the function V{t, x, y, z) differently in each 
subfamily. We have 

(r^ , V) with spherical symmetry; 
{u,w) = ^ {z,V) with plane symmetry; 

{x/y,V/y) with hyperbolic symmetry. 



2 dof 2 I 2 I 2 

r = X + y + z . 



(2.3) 



The formulae for matter density and pressure are known 
for each case, but will not be used in the present paper, 
so they are not quoted; see Ref. @. 

The Weyl tensor is proportional to /(u), and so with 
f(u) = the models given by (|2.2p become conformally 
flat. Then they become subcases of the Stephani solution 
given below, but are still more general than FLRW, see 
Ref. §\. 

The conformally flat Stephani solution |^ has the 
metric given by (|2.ip . the coordinates are still comoving, 
but the function V{t,x,y, z) is given explicitly by 



V 



^l + h{t) [ix~xo{t)f 



+ {y ^ yoit)f + (Z ~ Zoit)f]} , (2.4) 



where {R,k,xo,yo,zo) are arbitrary functions of t. This 
is easily seen to be a generalisation of the whole FLRW 
class, to which it reduces when (k,Xo,yo, zq) are all con- 
stant. The constants {xQ,yo,zo) can then be set to zero 
by a coordinate transformation, the constant k is the 
FLRW curvature index, and R{t) is the FLRW scale fac- 
tor. In general, the solution ()2.4p has no symmetry. 

As with the type D models, the formulae for mass den- 
sity and pressure are known, but will not be used here, 
see [1,11]. The mass density depends only on t. 

The parametrisation used in (|2.4p follows the original 
source [g and is designed so that the most popular rep- 
resentation of the FLRW models is easily obtained from 
it as a spatially homogeneous limit. It suggests that V is 
either quadratic in (a;, y, z) or does not depend on them. 
In fact, the metric (|2.ip is still conformally flat with 



V = Ao+ Agix"^ +y'^ + z"^) + Aix + Aay + A3Z, where 
the A's are arbitrary functions of t. However, even when 
Aq = 0, the quadratic terms can be restored with use 
of the Haantjes transformation, see Refs. 1,113 for the 
definition and examples of transformations. 

As already mentioned, all these solutions have zero 
shear and zero rotation. The quantity that makes them 
more general than FLRW is acceleration, which is in ev- 
ery case proportional to the spatial gradient of the goo 
component in the metric (j2.ip . Thus, the invariant con- 
dition for each of the metrics discussed above to reduce 
to an FLRW limit is 



d 

dx^ 



^1=0, 1 = 1,2,3. 



(2.5) 



Note that the metrics (|2.1I) do not allow a static limit 
- y,t must be nonzero. Static solutions that are shear- 
free and rotation- free obviously exist, but they form a 
separate branch in the family of solutions of Einstein's 
equations and cannot be recovered from (j2.1[) . 

There exists a very large body of literature on the 
SFN models, their various subcases and generalisations 

, but it is almost exclusively devoted to finding exam- 
ples of solutions of (12. 2p and of its charged generalisation, 
w,uu ~ ,fiu)w'^+g{u)w^. A general solution is not known. 



Part I 

The spherically symmetric models 

III. THE GEODESIC EQUATIONS 

For the spherically symmetric SFN models we first 
transform the spatial coordinates of (|2.1I) to the stan- 
dard spherical coordinates. Then the geodesic equations 
in the affine parametrisation are: 



dH 

ds2 



Eli 

F 



V 



1 



FWV^t 

' V,tu 



drV /di9 
ds] 



u sin d — — 



dip 
dJ 



V 



— + 2rFV„(y,,, 
_2^ — — -r 

y dsds 



dt du 
ds ds 

V,t v,u 



0, 



(3.1) 



V 



d§ 



dip 



+2r 



V 



dd 



dtp 
dl 



= 0, 



(3.2) 



d'^d ^V,tdtd^ 
ds^ V ds ds 



u y / ds ds 



3 



(3.3) 



(3.4) 



Equations p.ip - p. 21) can be simplified when we use 
the null condition: 



ds y ds ds \u V J ds ds 

+2coti9- — f- = 0. 
ds ds 



drV fddV ■^J'^vV ,^.r.2fdt 



What remains of them is: 

dH , fF.,t V,t , V,tt\ fdt 



(3.5) 



F V V,t \d 



^ V,t V dsds~ ^' 



(3.6) 



_ + ^(2,y,„-y„) _j -2--- 



r V y Mds 



0. 



(3.7) 



From p.3p - p. 41) it follows that the null geodesies 
are plane, i.e. that the coordinates can be adapted to 
each single geodesic so that it has = 7r/2 all along. 
However, in the following we will consider a bundle of 
geodesies emanating from a common source, and such an 
adaptation of coordinates would not be useful for that. 

In this and the next two sections, the same method as 
in Ref. yy] will be used, so, for readers' convenience, the 
relevant short excerpts from there are copied here with 
suitable modifications. 

For further calculations it is more convenient to use 
the coordinate r as a parameter, which will be non-affine. 
This is allowed, but with some caution. It is easily seen 
from p.7p that a curve obeying it, on which dr/ds = 
over some open range of s has dt/ds = in that range, 
and so is spacelike (the second possibility, 2uV,tu — V,t — 
0, could be compatible with (|2.2p only in the prohibited 
case V,t= 0). However, p.6p - p.7p do not guarantee 
that d7-/ds 7^ at all points; isolated points along a null 
geodesic, at which dr/ds — can exist. Thus, r can 
be used as a parameter on null geodesies only on such 
segments where ds/dr > or ds/dr < throughout. 

We have, for any coordinate: 



d^x" d^r dx" 
dr^ ds^ dr 



(3.8) 



Then, from p.7p we have: 



ds2 



{2uV,tu -V,t ) 



(3.9) 



Consequently, p.6p and p.3p - (|3.4I) become, using (13.91) 



dr^ 



d^ 
dr^ 



dV 
dr^ 



4r- 



it 



l\ dt 
dr 



-^(2«F„„-F„)(^ 
r V dr 



F Y,t 

3 

= 0, 



(3.10) 



d-d 
dr 



F'^V.t _ . ( dtY 1 



cos ■& sin -d 



dlyj 

dr 



-0, 



(3.11) 



dip 
dr 



2cot?9 



!^(2.y„„-y„)f^ 

r V dr 



dd dip 
dr dr 



(3.12) 



IV. THE REDSHIFT EQUATIONS 



Consider, in the metric (|2.ip . two light signals, the 
second one following the first after a short time-interval 
T, both emitted by the same source and arriving at the 
same observer. The equation of the trajectory of the first 
signal is 

(t,i^,^) = (T(r),e(r),<i>(r)), (4.1) 

the corresponding equation for the second signal is 

(i, z9, ^) = (r(r) + r(r), e(r) + C(r), $(r) + ^(r)). (4.2) 

Thus, we have to allow that while the first ray inter- 
sects the hypersurface of a given constant value of the 
r-coordinate at the point (i,-!?, (^) = (T, 9,$), the sec- 
ond ray intersects the same hypersurface at the point 
(i, 1?, 1^9) = (T -f T, -I- $ -t- V'). In general, it arrives at 
this hypersurface not only later, but also at a different 
spatial location. Thus, those two rays will not intersect 
the same succession of intermediate matter worldlines on 
the way. Note that, since the coordinates used here are 
comoving, both the source of light and the observer keep 
their spatial coordinates unchanged throughout history. 
Given this, and given that a pair of rays emitted by the 
same source and received by the same observer is con- 
sidered, (CV") = (0;0) holds at the point of emission 
and at the point of reception. However, the second ray 
is in general emitted in a different direction than the 
first one, and is received from a different direction by 
the observer. This means that in a general cosmological 
model the observed objects should drift across the sky. 
(See a brief quantitative discussion of this effect in the 
Lemaitre - Tolman model in Ref. [l|.) The directions 
of the two rays wiU be determined by (d9/dr, d<i>/dr) 
and (d0/dr + ^(r), d$/dr -|- 77(r)), respectively, where 
^ — dC/dr, r\ = dtp /dr. It will be assumed here that 
(dr/dr, C, ij), ^, rf) are small of the same order as t, so all 
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terms nonlinear in any of them and terms involving their 
products will be neglected. 

In writing out the equations of propagation of redshift, 
the symbol A will be used. It denotes the difference 
between the relevant expression taken at + r, r, i9 + 
(^,ip + ip) and at (t, r, -d, ip), linearized in (r, (, i.e. the 
difference between the value of a given quantity along the 
second ray and along the first ray, taken at a hypersurface 
of a given value of the parameter r (and, automatically, 
given u). For example F{t + T) - F{t) = /\F + ©(r^), 
V{t + r,u) - V{t,u) = Ay + ©(r^), A(di?/dr) = C 
This operation is a generalisation of the calculation by 
which Bondi fll| derived the redshift equation for radial 
null geodesies in the Lemaitre - Tolman model, see an 
account of that method in Ref. 111. 



Applying the A operation to p.lOp - p.l2p we obtain 
the following general equations of redshift propagation in 
the spherically symmetric SFN models:^ 



dr^ 



V,t 



dr 



F^tt F,t'^ V,ttt V,tt^ 
F 



[(2F,t V,t +FV,u ) {'^uV.tu -y.t ) 



F^ V^\^Ar_ 
V,t r dr \F V,tJ dr dr 



3F^V, 



r \ dr / dr 



(4.3) 



dr 



F^V, 



^(2.K,„-F,.) I +1 



[{2F,t V,t +FV,tt ) {2uV,tu -V,t ) 

2 

+ FV,t{2uV,ttu-V,tr '^'^ 

'L2l^2uV,tu-V,t)- — — 

r dr dr dr 

-cos(2i?)f^^ C-sin(2i?)^77 = 0, 
V dr / dr 



dr / dr 
dt d-0 dr 



(4.4) 



dr] 
dr 



-^{2uV,tu-V,t) 
r V dr 



--[{2F,tV,t+FV,u)i2uV,tu~V,t', 
r 



+ FV,t {2uV,uu 
2F^V,t 



{2uV,tu - 



' dr dr dr 



c 



„d?9 
-2cotz?— 
dr 



2cot?9^^ 
dr 



0. 



2 dddip 
sin^ d dr dr 

(4.5) 



Along with these, we can use also the result of A act- 
ing on the null condition (|3.5p . transformed to the r- 
parametrisation, it is 



(2FF,t V,t' + 2F^V,t V,t 



2A + ,sin(2^) 
dr \ dr 



■2F^V/ 
dip 



dt dr 
dr dr 



C + 2'usin^i9-^?7. (4.6) 
dr 



If T is the period of the electromagnetic wave measured 
in the rest frame of the source or of the observer, the 
redshift is given by 



T(robs) 
T(rom) 



1 + z(»'cm), 



(4.7) 



where the labels "obs" and "em" refer to the events of 
observation and emission of the ray, respectively P, [l3] . 
The period in the rest frame of an object is measured in 
the units of proper time of that object. In the metrics 
(12. ip the differential of the proper time, ds, is related 
to the differential of the coordinate time, dt, by ds = 
y/goodt, where 500 = {FV,t /V)'^- Thus, taking the T(r) 
found from (|4.3p - (|4.6I) as corresponding to the period 
T, we calculate the redshift from 



1 + z(rom) 



tobs.robs 



(4.8) 



Note that eq. (|4.3I) is decoupled from the other two and 
determines the redshift independently of (C, ip, ^, r]). This 
means that even for a non-radial ray the redshift changes 
with r by the same law as for a radial one. (Indeed, 
p.6p is the same for radial and non-radial rays, so A 
acting on it gives the same result in both cases.) This 
is a consequence of spherical symmetry and zero shear 
in the SFN models. Using eqs. (4.1) - (4.5) and (5.3) 
- (5.12) of Ref. Q one can verify that in the Szekeres 
model the equation (5.12) that determines r does depend 
on (CiV'iCi'?) even in the spherically symmetric subcase, 
where E,r = and shear is nonzero. 

Since C, = ip — Q at the observer, these quantities are 
not in fact observable. However, ^ and 77 are in general 
nonzero at the observer, which implies the change of di- 
rection toward the source with time. 



^ A quick way to calculate 1 14.311 - 114.61 1 is to take the differen- 
tial of the corresponding quantity at constant r and replace 
(dt,di9,d</3,d{di?/dr),d(d>^/dr)) by (r, C, ■0, f?)- 



V. REPEATABLE LIGHT PATHS (RLPS) 

We say that the light paths are repeatable when the 
rays sent between a given source and a given observer at 



5 



different times always proceed through the same sequence 
of intermediate particles of the cosmic medium. This 
means, when the rays are registered at an r-hypersurface 
of coordinate radius r, they arrive there at the same spa- 
tial location (only at different time instants), i.e.: 



C = = ^ = - 



(5.1) 



all along each ray. Substituting (|5.ip in (|4.4p - (|4.6I) one 
obtains the conditions that have to be obeyed in order 
that RLPs exist. We first define: 

Xit, u) [{2F,t V,t +FV,u ) {2uV,tu -V.t ) 

+ FV,t {2uV,ttu -V,tt )] + 2FV,t {2uV,tu ~V,t ) ^ 

or 



and then the conditions are 
F dtdz? 

— ^;r:r = 0' 

r ar ar 

F dtdip „ 

— ^A~ir ^ ^' 

r ar ar 



2FV,t 



{F,t V,t +FV,u ) + FV,t ^ 
ar ar 



dr 



(5.2) 

(5.3) 
(5.4) 
0. (5.5) 



Discarding in (j5.5p the impossible solutions dt/dr = 
(this would be a spacelike curve) and FV,t — (prohib- 
ited - see (|2.ip ) we obtain: 



dr 
dr 



F 



V,t 



di^ 
dr 



(5.6) 



Equations (|5.3p - (|5.4p imply that either d-d/dr = 
d(p/dr = or X = 0. The first case defines a radial null 
geodesic. This means that radial null geodesies are RLPs, 
as expected in a spherically symmetric model. Then (15. 6p 
just defines the redshift as a function of r. 

To check the other possibility we substitute (|5.6p in 
X ~ 0. The result is, after discarding the factor 
2u (F^/r) (di/dr)T: 



V,tt V,tu -V,t V,uu = 0. 
A general solution of this equation is 

V^a{u)S{t)+/3{u), 



(5.7) 



(5.8) 



where a, /3 and S are arbitrary functions of their respec- 
tive arguments. This must be compatible with (|2.2p . We 
disregard the cases aS,t — because they lead to the 
impossible condition V,t = - see the penultimate para- 
graph of Sec. HB Then, (jO)) is compatible with ([^ 
only when / — 0, i.e. when the metric is conformally 
flat. The conformally flat case belongs to the Stephani 
class (|2.4p , which will be discussed in part IV of this pa- 
per. However, in that class the spherically symmetric 



case emerges among many others, and it will be more 
convenient to discuss it here.^ 

When / = 0, (|5.8p substituted in (|2.2p gives = 
/3,„n = 0, i.e. 



V^AiS + Bi + {A2S + B2)r^. 
This is in the Stephani class (|2.4p . with 



(5.9) 



1 



R{t) 
xo = yo 



= Zo 



= 0. 



A2S+B2 

AiS + Bi' 



(5.10) 



An FLRW limit results from (|5.9p when the zero- 
acceleration condition (|2.5p is fulfilled, i.e. when A2B1 — 
A1B2 = 0. If A2 = 0, then in the FLRW limit A1B2 = 0. 
But ^1 = ^2 = is the prohibited case V^t = 0, while 
^2 = ^2 = is the fc = FLRW model. So, for a general 
FLRW limit of dO]) we have A2 7^ and 



Bi = A1B2/A2. 



Then, in the FLRW limit 



V 



m 



1 -I- -r' 



(5.11) 



(5.12) 



where 



R 



Ai {A2S + B2) 



k,'-^. (5.13) 



Since S is arbitrary and (^1,^2) can have any signs, 
this shows that the model defined by (|5.9p reproduces 
the whole FLRW class in the zero-acceleration limit. 

In the solution given by (|5.9p - (I5.10p all null geodesies 
are RLPs, since the metric obeys (|5.7p with no conditions 
imposed on the vectors tangent to null geodesies. 

Thus, the result is: 

Corollary 1: In a general spherically symmetric 
SFN model the only repeatable light paths are radial 
null geodesies. In the subcase defined by (|5.9I) all null 
geodesies are RLPs. This subcase contains the whole 
FLRW class, but is more general than FLRW because it 
has nonzero acceleration. It has zero Weyl tensor, but is 
less general than the spherically symmetric limit of the 



^ It may be verified that 1 15.61 1 is indeed a first integral of 114.31 1 
modulo the RLP conditions. Hint: take a derivative of 115.61 1 
along a null geodesic (defined by TD/dr = (dt/dr)d/dt + d/dr), 
then use 113.101 1 and 115.61 1 to eliminate d-^t/dr^ and dr/dr. The 
result will be equal to 114.31 1 provided that: 



dt 



2 (V,tt V,tu -V,t V,ttu ) 



FV^Y 
dr J 



: 0. 



The expression in ( ) vanishes when (15.71 1 is fulfilled, the ex- 
pression in [ ] vanishes in virtue of the null condition when the 
geodesic is radial. Thus, l|5.6| l is a first integral of ||4.3| I when any 
of the RLP conditions is fulfilled. 
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Stephani solution (j2.4p : in the latter k{t) and R{t) are 
two independent functions. 



Part II 

The plane symmetric models 

VI. THE GEODESIC EQUATIONS 

The scheme of the calculation is here the same as in 
Sections IIIII -IVl only the explicit forms of the equations 
are different. Thus, we will limit the explanation and the 
presentation of intermediate expressions to a necessary 
minimum. 

Analogously to the spherical case we note that z 
can be used as a parameter on open intervals of each 
geodesic and use the d^z/ds^ geodesic equation to carry 
out the reparametrisation.'^ The geodesic equations 
parametrised by z are then: 



V,tz dt 



V,t dz 



d^x 



dz"^ \ dz 



dt \ " da; 



dz 



= 0, 



dz^ ' ' \ dz / dz 



(6.1) 
(6.2) 
(6.3) 



Equations (|6.2I) and (|6.3p show that (ydx/dz—xdy/dz) is 
a constant of the motion. This means that the projection 
of each geodesic on the (z, y) plane is a straight line. 
In obtaining the redshift equations via the A opera- 



tion defined in Sec. IIVI we take (C, V') '== ^ix,y) and 



(^,77) =^ (d/dz)(^, -0). Since we will not use the redshift 
equations in full, we do not display them and proceed to 
the conditions for the RLPs. 



VII. REPEATABLE LIGHT PATHS 

We substitute ^ = 'i/' = ^ = '/ = Oin the redshift equa- 
tions and in the null condition, and obtain the conditions 
for RLPs: 

X{t, z) = [(2F,t V,t +FV,u ) V.t. +FV,t V,u. ] 



^ If dz/ds = on an open interval of the geodesic, then either (a) 
dt/ds = 0, which can hold only on a spacelike geodesic, or (b) 
V,t = 0, which is a prohibited case (see the remark below 1 12. 5t l. 
or (c) V,tz = 0, which is compatible with I I2.2I I only when V,t = 
or / = 0. The latter case is the conformally flat model that will 
be dealt with in Sec. IVIII 



dr 

+2FV,tV,t.—, 

^ dt dx Q 
dz dz 

dz dz 



2FV,t 



dt 



dt 



iF,t V,t +FV,u ) - 



FV^t 



dr 



(7.1) 
(7.2) 
(7.3) 
0.(7.4) 



Discarding the impossible cases, there are two sets of 
solutions to ([7^ - ((73)) : (a) dx/dz = dy/dz = - the 
analogues of radial null geodesies; (b) x = - the subset 
of the plane symmetric models in which all null geodesies 
are RLPs. Using (17. 4p in x = to eliminate dr/dz, we 
obtain: 



V,u V,t. -V,t V,ttz = 0, 
a general solution of which is 

V = a{z)S{t)+f3{z), 



(7.5) 



(7.6) 



where a, /3 and S are arbitrary functions of their respec- 
tive arguments. Just as (|5.8p . this is compatible with 
(|2.2p only when the metric becomes conformally flat. 
Then a,zz = 13, zz = and 



V ^AiS + Bi + {A2S + B2)z. 



(7.7) 



This case is not covered by the parametrisation of (|2.4p . 
but can be brought into the form of (|2.4p by the following 
special Haantjes transformation; 



{x,y,z) 



(x',y',z' + c(x'Vy'' + z'2)) 



T 



r = l + 2Cz' + C^(x' 



2 , a , /2 
y + z 



where 
(7.8) 



where C is an arbitrary constant, a group parameter of 
the transformation (see Refs. 0, for hints on how to 
handle the calculations). The result of 



V = C[C {AiS + Bi) + A2S + B2] (x'^ + ?y'' + z'' 

+ [2C {AiS + Bi) + A2S + B2] z' + AiS + Bi. (7.9) 

For later reference we will need (|7.9p cast explicitly in 
the form (12. 4p . The functions appearing in (|2.4p are in 
our present case: 



1 

'r 

7 
k 



dct 



{A2S + B2y 

47 

{A^S + Bi) + C {A2S + B2) 
I672 



{A2S + B2y 



C {AiS + Bi) + -{A2S + B2 



(7.10) 

(7.11) 
(7.12) 

(7.13) 



Equation (|7.12p suggests that fc < necessarily. How- 
ever, the fc = case follows directly from (17.71) . it is 
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A2 = B2 = 0, and then 1/R = AiS + Bi. Thus, the 
plane symmetric subcase of the Stephani model can re- 
produce the fc < FLRW limits. 

The general prescription for the FLRW limit can be 
found by substituting ()7.7p in the zero-acceleration con- 
dition (|2.4p : the result is 



A1B2 ~ A2B1 = 0. 



(7.14) 



We will need the formulae for R and k parametrised 
by Zq rather than by S. From (|7.13l) we find 



S 



BiC + B2/2 - {BiC + B2) Czo 



{CAi + A2)Czo-CAi 
and then from (frT2l) and (TTTO)) 



A2/2' 



1 

2K 



dcf k 
^ = 4i?' 



(CAi + A2) Czo - CAi - A2/2 



C2 {A2B1 - A1B2) 



1 



2KR 



1 2 , fo 1_ 

2""" C 2C2' 



(7.15) 

, (7.16) 
(7.17) 
(7.18) 



Equations (|7.16|) and (|7.18p will be needed to recognise 
the plane symmetric subcase among the multitude of 
cases discussed in Appendix A. 

Note that the form (|7.16p - (|7.18p does not allow taking 
the FLRW limit (|7.14p . This is logical, since in this limit 
Zq becomes constant while K and R do not, so zq cannot 
be used as a parameter. 

We thus have: 

Corollary 2: In a general plane symmetric SFN 
model represented as in p.l|) - p.3|) the only repeat- 
able light paths are null geodesies on which x and y are 
constant. In the subcase defined by (|7.7p . which includes 
the k < FLRW subclass, all null geodesies are RLPs. 
This subcase is conformally flat, but less general than the 
conformally flat limit of the plane symmetric case. 

The condition fc < is consistent with what is known 
about the relation between plane symmetric models and 
their FLRW limits 

Part III 

The hyperbolically symmetric 
models 

VIII. THE GEODESIC EQUATIONS 

There is a certain complication in discussing this case: 
either the metric can be similar to (j2.1l) or the equation 
similar to (|2.2p . but not both things at once. We choose 
to make the metric similar. We use eq. (2.1) of Ref. 

0, rename the functions by {30, Y, b) '= {1/F, V, /) and 
transform the r-coordinate used there by r = exp(r'). 



Dropping the prime, the metric then becomes: 

2 



ds^ ^i^] dt^ - -L (dr^ + d^2 + sijV ^^^2) ^ 



V 



and the function V{t, r) must obey: 



(8.1) 



(8.2) 



See [i| for the transformation from to - ([23]) . 

Just as in Sec. IVII we will now skip most of the ex- 
planation and of the intermediate expressions because 
the calculations exactly parallel those for the spherically 
symmetric models presented in Sections IIIII - |Vl 

Using the (Pr/ds^ geodesic equation we change the 
parameter to the non-affine r and obtain: 



dr^ 











V,t) 


(dr) 



F'V.,V.,A%) =0, 



dr^ dr y dr / K'^'"' J 



(8.3) 

2 



= 0, 



(8.4) 



- F'V,tV,tr-^[-r] +2cothi?— =0. 
dr^ dr \ dr / dr dr 

(8.5) 

We skip through the general redshift equations 
and proceed to the RLP equations, where this time 

(C, V) =' A(^,¥') and (e,77) (d/dz)(C, V)- 



IX. REPEATABLE LIGHT PATHS 



When C = V' = C 
become: 



0, the redshift equations 



X{t, r) [(2F,t V,t +FV,u ) V,tr +FV,t V,ur ] 



-2FV,tV,tr'^. 

dr 



dr dr ' 
dr dr 



2FV,t 



As before, (|9.4I) implies 

dr _ 
d^ ^ " 



dt 



{F,t V,t +FV,u ) -^T + FV, 
dr 



dr 
dr 



di 
dr 



F,t , V,t 



F V,t J dr 



dt 



(9.1) 
(9.2) 
(9.3) 
0. 

(9.4) 



(9.5) 
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Equations (H^) - tell us that either d-&/dr = 

d(p/dr — or X = ^- The first case defines a pseudo- 
radial null geodesic. The other possibility is x = 0. Using 
(1^3)1 we get 

V,uV,tr-V,tV,Ur=0. (9.6) 

A general solution of this equation is 

V ^ a{r)S{t) + /3{r), (9.7) 

where a, /3 and S are arbitrary functions of their respec- 
tive arguments. As before, this is compatible with ()8.2|) 
only when f — and the solution becomes conformally 
flat. Then, from ((8?2|) . V must obey, V,rr+V = 0, and 
the solution is 



V = {AiS + Bi) sinr + {A2S + B2) cosr. 



(9.8) 



This is transformed to the standard parametrisation (|2.4p 
by the following chain of transformations. 

Wc first transform the 2-dimensional metric {dd^ + 
sinh^ -ddip^) into (dr^ -l-e^'^dz^). An explicit prescription 



for this is given in Appendix A to Ref. [9| . This does not 
affect V because d and tp are not present in it. Then we 
carry out the transformation 



arcsm ■ 



y 



The metric then becomes: 



r = --ln(a;2+y2). (9.9) 



ds2 



fFW,t 



\ W 



de 



dx^ + dy"^ + dz^ 



(9.10) 



where 



W = .Jx^+yW = [AiS + Bi)y + {A2S + B2) x. 

(9.11) 

The metric (|9.10p - (19.111) is within the Stephani class, 
as can be seen by carrying out the following Haantjes 
transformation: 

_ {x' + C(x'^ +y'^ + z'^),y',z') 
^' = l + 2Cx' + C^{x''+y'^ + z'') ■ ^^-^^^ 

After this, the metric acquires the form (|2.ip with V 
replaced by 

Wi = C {A2S + B2) (x'^ + y'^ + z'^^ 

+ {A2S + B2)x' + {AiS + Bi)y'. (9.13) 
When Wi is cast in the form of (|2.4I) . we obtain 



2/0 
K 



_ 1 

^ ^2' 

lAiS + Bi 

~ ^2A2S + B2' 

dof fc 



4i? 



C{A2S + B2)., 



(9.14) 
(9.15) 
(9.16) 



1 _ C {^A^s + B^f + {A2S + B2f 

4 ^2^ + 52 ■ ^ ^ 

The constant a;o can be transformed to xq = 0. 

The FLRW hmit follows when A1B2-A2B1 = 0. From 
(|9T6)1 - (I9T71) we find 



16{A2S + B2Y 



(9.18) 



(Ai5 + Si)' + (A2^-f S2)" 



which shows that fc < necessarily. The V given by ^ 
cannot be made independent of r, so the case A: = is 
not contained in this model. 

For later reference we will need the formulae for K and 
R parametrised by j/o rather than by S. We have 



S = - 

1 _ 
2K " 
1 

2KR 



Bi + 2B2yo 
Ai + 2A2yo ' 

Ai + 2A2yo 
2C (A1B2- A2B1)' 
1 o 1 



2C 



Vo 



8C' 



(9.19) 
(9.20) 
(9.21) 



Similarly to the plane symmetric model, in this form the 
FLRW limit A1B2 — A2B1 — cannot be taken, for the 
same reason: in the FLRW models yo is constant, while 
K and R are not, so yo cannot be used as a parameter. 
Thus we have: 

Corollary 3: In a general hyperbolically symmetric 
SFN model the only repeatable light paths are the ana- 
logues of radial null geodesies, on which d and (p in (|8.ip 
are constant. In the subcase defined by (|9.10p - (|9.1ip . 
which includes the FLRW models with k < 0, all null 
geodesies are RLPs. This subcase is conformally flat, 
but less general than the conformally flat limit of (|8.ip . 

Part IV 

The Stephani model 

X. THE GEODESIC EQUATIONS 

The metric of this model is (j2.ip . with V given by (|2.4I) . 
We introduce the abbreviation: 



D 1^' Fit)V,t IV. 



(10.1) 



The geodesic equations in the affine parametrisation, 
with the null condition already incorporated, are: 

— , dx D.ydy D,^dz\ dt 

ds2 \ D ds D ds B ds) ds 



(10.2) 
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V ds 



V,y dy V,z dz 
V ds 



ds2 



ds2 



V ds 

2 



V^dt_ 
V ds 



df 
ds, 

V^dx 
V ds 

dt' 

V,x da; 
~dJ 



V^dt 

y ds 



dx 
d^ 

2 



= 0, (10.3) 



- 2 
+ i^"^ 

- 2(^^ + ^^ + ^^)^ = 0, (10.4) 

+ F^V,tV,t. 

- 2 



dy 
ds 

V^dz\ dy 
y ds y ds 

2 



dy 



dzV 
dlj 



y ds 



dz 
dl 



= 0. (10.5) 



Note that if any of dx'/ds, i = 1,2,3 (where (x^,x2,a;'^) 

'= {x, y, z)) is zero along an open interval of a null 
geodesic, then piH)) - (ITU3)) imply either dt/ds = 
in the same interval (which is impossible on a null curve) 
or V,ti = for the respective i. This second possibility 
must be investigated. Suppose first that 



V,t. = V,ty = V,t. = 0. 



This means: 



2^ ) ,t = a;o,t = yo,t = zo,t 0. 



(10.6) 



(10.7) 



The metric is then spherically symmetric and the con- 
stants {xQ,yQ,zo) can be set equal to zero by coordinate 
transformations. This case was dealt with at the end of 
Sec. El so we disregard it here. 

Hence, we may assume that at least one of the quanti- 
ties in (110. 6p is nonzero. Since a general Stephani metric 
does not change its form under any permutation of the 
{x, y, z) coordinates (accompanied by a suitable renam- 
ing of the (xo, ya, zq) functions), we may assume without 
loss of generality that V,tz ^ 0, and then (jlO.Sp shows 
that dz I ds cannot be zero over an open interval of a null 
geodesic. Consequently, z can be chosen as a (non-affine) 
parameter, with an analogous cautionary remark to the 
one given above (13.81) . 

Consequently, using (|10.5p and (|3.8p for d^z/ds^ we 
change the parameter to z in (|10.2p - p0.4p and obtain: 



d^i 
d^ 



ax 
d^ 

dz2 



dz) '^{f ^ V,tJ \dz 

(10.9) 
(10.10) 



XI. REPEATABLE LIGHT PATHS 

In this case we skip the redshift equations because they 
are complicated and voluminous, while we are not going 
to make any direct use of them. In this section, the 
meaning of the symbols in the A operation is 



A(t, X, y , dx/ dz, dy/dz) =^ (r, C,fp,t v) 



(11.1) 



We proceed to the RLP conditions that are obtained by 
applying the A operation to ()10.9p - (jlO.lOp and immedi- 
ately assuming ( = ip — ^ — ri = 0.^ We do not consider 
the result of A acting on (jlO.Sp because it is the equation 
for T that will define the redshift propagation along the 
emergent RLP, and it does not lead to any limitation on 
the metric. 

The null condition in the z-parametrization is: 



dt 



dz 



da; 



dz 



and the RLP condition resulting from it is: 
, dt dr 

{F,t v,t +FV,u + :r = 

dz dz 



(11.2) 



(11.3) 



(we ignore the cases FV,t = - prohibited in (|2.ip . and 
dt/dz = ~ which defines a spacelike curve). 
We denote: 



„ dcf V,ti 



„ dcf „ 



1,2,3, (11.4) 



where (a;-'^, a;^, a;'^) ^= (a;, y, z). We find dr/dz from (jll.Sp 
and use it in the RLP conditions resulting from (|10.9p - 
(jlO.lOp in the way described above. The result is: 



da; 



n dy 



0, 



0. 



(11.5) 



(11.6) 



Since the analysis of (jll.Sp ~ (|11.6p is complicated, we 
state here only the results, and present the calculations 
in Appendix S The RLPs defined by (fTT31) - pT^ 
exist in the following cases: 

(1) (Case 1.1.1.1.2.1.1 of AppendixEJ 

When yo is defined by (|A23]) . K =^ /c/(4i?) by (|A24ll 
and R by (|A39p . some of the null geodesies are RLPs, 
and they are solutions of (jA52p . with T\ and Q\ given 
by (|X40)) - (|M2)) and ((X50)) - ((X5T|) . It is not known 
whether this subcase of the Stephani solution admits an 



^ Similarly to what footnote #1 says, a quick way to obtain the 
equations that follow is to take the differentials of l|10.8|l - Ill0.10| l 
at constant (x, j/, z) and replace At by r. 
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FLRW limit. If it does, then only after a reparametrisa- 
tion. As noted for the plane and hyperbolically symmet- 
ric cases, when k and R are parametrised by xq, uq or zq, 
the FLRW limit cannot be calculated. 

Several subcases of this spacetime appear separately in 
Appendix \^ but they are not listed here. 

(2) (Case 1.1.1.2.1.2.1 of AppendixEJ 

When Uq — Dixq, zq — C3X0, and R is determined by 
(|A64p . but K{t) = k/{AR) is arbitrary, again some of the 
null geodesies are RLPs. They are determined by (jA52p . 
where Fx and Qi are given by (jA66l) ~ (|A67p . 

Also here, several subcases appear separately, but they 
do not admit more RLPs. 

(3) (Case 1.1.2.2.2 of AppcndixEl 

This is a copy of the subcase of the spherically sym- 
metric Stephani solution discussed in Sec. |Vl All of its 
null geodesies are RLPs, and it contains the FLRW limit 
in full generality. It is not identical to FLRW because in 
general it has nonzero acceleration. 

(4) (Case 1.2.1.2 of AppendixEJ 

Then xo = CiZq, j/o = Ca-^O: R{t) is arbitrary, V and k 
are determined by (jAllip and (jA112[) . This is an axially 
symmetric subcase of the Stephani model, and its RLPs 
are determined by (jA113p . All the RLPs intersect the 
symmetry axis a; = y = in each space of constant t. 

(5) (Case 2.1 of Appendix El 

Then j/q — C2X0, zq — C3X0, K — k/{AR) is deter- 
mined by (IA137I) and R is determined by (|A141|) . As 
indicated there, C2 = C3 = may be achieved by a co- 
ordinate transformation, and then the model is seen to 
be axially symmetric. All of its null geodesies are RLPs. 
It is less general than the previous one (see under (4) 
above) because the former has R{t) arbitrary, while the 
current one has R determined by (jA141l) . 



corresponding conformally flat limits of the SS, PS and 
HS cases. 

For the Stephani models the situation is summarised 
in Sec. IXII In general, no RLPs exist. The conditions 
of existence of RLPs, (lll.5|) - (|11.6p . put limitations on 
the Stephani metric. There exist subcases in which some 
of the null geodesies are RLPs, for example the axially 
symmetric subcase given by (|A111I) - (jA112p . There ex- 
ists also a subclass in which all null geodesies are RLPs, 
it is axially symmetric as well. 

This study, as explained in the introduction, was moti- 
vated by the question: what is the geometrical condition 
for the existence of RLPs; is it vanishing shear in the flow 
of the cosmic medium, as suggested by the result of Ref. 
[l[ , or rather a high symmetry of the spacetime? The so- 
lutions of Einstein's equations investigated in the present 
paper all have zero shear, and yet their null geodesies 
are RLPs only in special situations. Thus, the key to 
the RLP property is a symmetry of the spacetime rather 
than vanishing shear. 



Appendix A: RLPs obeying (|11.5p - (|11.6p 

The calculations in this Appendix are trivial in prin- 
ciple. The reason why they are presented in some detail 
is that the various separate subcases form a very com- 
plicated binary tree that would be difficult to duplicate 
without a guidebook. 



Case 1: The general case: ^ 
Then (llTSl) - (fTTB)) become 

dx _ Ay _ Gy 

dz ~ G ' dz~ G' 



(Al) 



XII. SUMMARY 

The existence of repeatable light paths (RLPs) was in- 
vestigated for the cosmological models found by Barnes 
[3| and Stephani They are called shearfree nor- 

mal (SFN) because the perfect fluid source in the Ein- 
stein equations moves with zero shear and zero rotation. 
The Barnes models are either spherically symmetric (SS) 
or plane symmetric (PS) or hyperbolically symmetric 
(HS). In general, in each of these classes only those null 
geodesies are RLPs that are orthogonal to the symmetry 
orbits (the radial ones in the SS case). However, each one 
contains a subclass in which all null geodesies are RLPs. 
This subclass is in each case more general than FLRW. 
In the SS case, this special subclass has the metric (12. 4p 
with V given by (15. 9p and (|5.10p , and contains the whole 
FLRW family. In the PS case, the special subclass has 
V given by ([7^ - ((71^ and contains only those FLRW 
models for which the curvature index fc < 0. In the HS 
case, the special subclass is given by (|9.13p - (I9.17P and 
contains only the fc < FLRW models. All these special 
subclasses are conformally flat, but less general than the 



We calculate the derivatives of dx / dz and dy / dz along a 
null geodesic by the rule 

D dt d dx d dy d d 
dz dz dt dz dx dz dy dz 



and substitute the resuhs in ([TII^ - (|10.10p . Then 
we use (jAip and (|11.2p to eliminate dx/dz, dy/dz and 
{dt/dz)'^. The resulting equations are 

Gz {GzGx,t ~ GxGz,t) — h Gx (GzGx,x ~ GxGz.x) 
dz 

+ Gy {GzGx^y ~ GxGz^y) + Gz {GzGx^z ^ GxGz^z) 

+ (HxGz - HzGx) {Gx" + G/ + G,2) ^ 0, (A2) 

Gz {GzGy^t ~ GyGz,t) + Gx {GzGy^x ~ GyGz,x) 

+ Gy {GzGy,y ~ G yG z ,y) + Gz {G z.G y , z " G yG z . z) 

+ [HyGz - HzGy) {Gx^ + Gy" + G, ^ ) - 0. (A3) 



After (|11.4p are substituted for Gi and the terms free 
of dtjdz sum up to zero in each of (jA2p - (jA3l) . so both 
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coefficients of dt/dz must be zero, too. This implies: 
Gx \ ( G 



0. 



(A4) 



G/ \ 

The integrals of these are: 

Gx= Tx{x,y,z)G^_, Gy = Qx{x,y,z)G^, (A5) 

where Ti and Qi^ i = 1,2,3 are arbitrary functions. Re- 
caUing (jll.4p . these are integrated again with the result: 



Hx = TxH^ + Tiix, y, z), 
Hy = QiH^ +g2{x,y,z). 



(A6) 



Finally, recalling the definitions of and Hy from 
(lll.4p . eqs. (IA6|) are integrated with the result: 



V,x^TiV„+T2V + T3{x,y,z), (A7) 

v,y^giV„+g2V + g3{x,y,z). (as) 

We introduce the following conventions: 

• Ci, Di, Ci,di^ Ei,i = 1,2,... will denote arbitrary 
constants, 

• -^ii Sj, j = 4, 5, . . . will denote arbitrary functions of 
spatial coordinates, not necessarily of all of them. 
Since the alternatives considered below will be in most 
cases mutually exclusive, we will re-use the same names 
of constants and functions with different meanings. 

In every case we will follow the same scheme of reason- 
ing. Our initial equation (IE) will be (IA7|) or (|A8[) . usu- 
ally multiplied by some factor (1/27^ or or ^2^^)- 
Then, the following operations will be executed on IE, 
each one followed by conclusions: 

1) Differentiate IE by y and t; 

2) Differentiate IE by y alone; 

3) Differentiate IE by x and t; 

4) Differentiate IE by x alone; 

5) Differentiate IE by z and t; 

6) Differentiate IE by z alone. 

At this stage, {J^i,gi),i — 1,2,3, K and R will be de- 
fined, and the corresponding Stcphani model admitting 
RLPs (if any) will be identified. The functions (J^i,t/i), 
substituted in (jAl|) . will define the RLPs. 

We will present in detail the whole 6-step procedure 
only for the first case considered. 

It will turn out in several places along the way that the 
function xo{t) is in fact constant. In those cases, we will 
assume that a;o = because this result can be achieved 
by the coordinate transformation x — x' -\- xq. The same 
is true for the pairs {yo,y) and {zo,z). 

Case 1.1: T2 

Then (jA7p can be written as: 

X — xq , , T^ 1 



1 r 
2 



{z - zo) — 



T2 



2KR 



(A9) 



{x ~ xof + {y- yof + {z - zof 



2KT2' 



where we have introduced the symbol: 



K = . 

AR 



(AlO) 



In writing (jA9l) we assumed fc 7^ because fc = is the 
spatially flat FLRW model, in which we know that all 
null geodesies are RLPs. Taking the second derivative of 
^9[) by y and t we obtain: 



1 



x^^t \y-]^v- ^o.t I ;^ ) +yo,t - I 2if j ' " ' 



1 



(All) 



Case 1.1.1: xo,t 7^ 



Then we divide (|A11|) by xo,t and differentiate the re- 
sult by t, obtaining: 



zo,t \ 








1 




_Xo.t 


(^) 



yo,t 

xo,t 



(A12) 



Case 1.1.1.1: (zo,t/xo,t) 7^ 



Then we divide (|A12|) by (zo,t/a;o,t) ,t and differentiate 
the result by t. We get: 



xo,t 



xo,t 



,t-xW( J),y = 0, (A13) 



where: 



xit) 



dcf 



1 



1 



_xo,t \2K 
Case 1.1.1.1.1: x{t) + 



xq± 



(AM) 



Then we get: 



a;o,t/ / \xQ,t 
— C\— constant. 



T2 



(A15) 



both expressions being constant because the first one de- 
pends only on t, while the second one depends only on 
(x, y, z). Integrating both equations we obtain: 



^1 n 

2K = 



yo + DsZo + Di, (A16) 
T3 = {Ciy + Ti{x,z))T2. (A17) 



In principle, we would have to consider the cases Ci 7^ 
and Ci = separately. However, they lead to the same 
result. When Ci ^ 0, (|XT6)l d ete rmine s 1/{2K). Then 
we substitute (|A16P and (|A17p in (IA12P and obtain 



zo,t 

Xo,t 



^2 



'As 



(A18) 
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Since in Case 1.1.1.1 the first factor is nonzero, we have: 
^1 = {Dsy + Mx, z)) T2. (A19) 
When we substitute (UT7| and (|JT9)) in (|IlT|) we obtain: 
1 



ry +Di 



0. 



(A20) 



Since xo,t 7^ in our current Case 1.1.1, we have: 



1 



— = -Diy + Tq{x,z). 
Using dMU, (|JT9l) and ((MtI) in (jA9l) we obtain: 



(A21) 



1 1 



[x ~ xof + 2/^ + y^ + (2 - zo)' 



2KR 2 
~T,^D,y. 



The derivative of this by y now gives ~D\x — D^z + y + 
D4, which is a clear contradiction. 

When Ci = 0, K remains undetermined, and instead 
(jAlsp determines yo = —Dixq — D3Z0. Equations (|A19P 
and ((MT|) still follow and (|X22]) results with D4 = 0, 
leading to the same contradiction. 

Case 1.1.1.1.1 thus turned out to be empty, and we go 
back to (|A13|) to consider: 

Case 1.1.1.1.2: x{t) = 



Then (|A13p implies [{yo,t/xo,t) ,t / {zo.t/xo^t) ,t ] ,t = 0, 
which is integrated with the result: 



yo = CiXq + C^zq 



(A23) 



(the additive constant was set to zero by a transformation 
y = y'+ constant), and x{t) = implies, via (|A14p 



^ = Dixo + D3Z0 + D4. 



(A24) 



We substitute (jMSl) - (jXM]) in (|JT2|) and obtain 



9u '^-^-^(Ji'^+^s 



= 0. (A25) 



Since we are still in Case 1.1.1.1, where (zo,t/a;o,t) ,t ^ 0, 
the above leads to 

^1 - + (Mx, z) - Cay) T2. (A26) 

Substituting this and (|M3)) - (|X24l) in (PIT|) we get 



0. (A27) 



In Case 1.1.1, where xo.t 7^ 0, the above implies 

DiT3^iCiy + T5ix,z))T2-l. (A28) 



Case 1.1.1.1.2.1: Di ^ 



Then (jA28l) determines ^"3. Using (jA26|) . (|A24l) and 
(|M3)) in dMI we get: 



1 r 

'2 



{x - xo f +y^ + (Cixo + CazoY + 



{z - ZqY 



•r2 



(A29) 



We differentiate this by y alone and obtain 



3 I ^ I vy -C's 



(A22) — = 1. 



(A30) 

Integrating this back and using (jA28l) for we obtain 
1 ^D,y^ + {CiD3~C3D,)yz + CiD4y 



T2 Dix + D3Z + D4 

Using ([Ms]) and (|X3T|) in (|X29)) we obtain 

1 



Te{x,z). 
(A31) 



xTg^x, z) = (z- zq) Tiix, z) + — — — + a;oJ'5(x, z) 

ZiX ri 



(X -~ Xof + {CiXq + Cszof + (z - Zq)'^ 



+ — {D3Z + D4) iMx, z) - Ffi{x, z)) . (A32) 

Taking the second derivative of this by x and t we obtain 
zo.t 



XQ.t 



-^4,,x — 1 + ^5,x = 0. 



(A33) 



In Case 1.1.1.1, where (zo.t/a^o.t) ,t 7^ 0, this implies 

Fi = Ti{z), Fz=x + Tj{z). (A34) 



Putting this in (jA32p and taking its derivative by x alone, 
then integrating back, we get 



D3Z + D4 



1 2 , xiD^z + Di) 
Fq{x,z) = -X +- 



(A35) 

Using this in (IA32[) . then taking the second derivative by 
z and t we get 



za,t 

XQ.t 



{J'i,z + 1) + J'l^z = 0. 



(A36) 



In the present case this means 

J-T = Ci = constant, = —z + C^. (A37) 

Feeding this information in (|A32p and using (jA24p we get 



J^s(z 
1 



^z + \z'-C,z-^C, = E„ (A38) 



R 



(DiXo + D3Z0 + D^) + dxo - C5Z0 
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xq^ + (Cixo + C3Z0) + za 



El, 



(A39) 



the two expressions being constant because they are 
equal, while the first one depends only on z and the sec- 
ond one only on t. Equations (|A38|) - (|A39|) define R{t) 
and J-s{z), and are the final solution of (|A7|) . 

Picking up the pieces, we obtain the following formula: 



^1 



dcf 1 



(A40) 



^3 (3 



V 



z^) + iCiD3-C3Di)xy 



-Dixz + (QDa + C^Di)x - CsD^y - D^z 
+C5D^~D3Ei, (A41) 



= ^Di{x 



2 2\ 

V -z ) 



D3XZ 

{CiDs - CMyz + Dix + CiDiV 
iC4D3 + C5Di)z + CiDi + DiEi. (A42) 



Now we will deal with (|A8p . still within Case 
1.1.1.1.2.1. 

Case 1.1.1.1.2.1.1: ^2 7^ 



Then (|A8|) is written as 
y - C\xo - C3Z0 



\z - Zq) — 
SJ2 



2KR 



j_g3 

2X02 

2 



(a; - xo)^ + {y - Cixo - Cszof + {z - zq) 



(A43) 



We follow exactly the same sequence of steps that we did 
in solving ()A9p . We take the second derivative of (|A43P 
by y and t and obtain, using (jA24l) : 



zo,t 
'xo,t 
zo,t 



ZQ,t 
XQ.f 



(A44) 



XQ.t 



zo,t 



Since we are in Case 1.1.1.1, the coefficients of (^o,t/2;o,t) 
and the remaining terms must balance separately. More- 
over, we are in Case 1.1.1.1.2.1, where Di ^ 0. Integrat- 
ing the two equations with respect to y we obtain 



Q2 

Q2 



+D^g^{x,z), 

Ci ( 1 



^ ) +^4(2;, z). 



iyA'-Q-2 



Gbix^z). 



(A45) 
(A46) 



We substitute ([HS]) - ([Xiej) in (jXIgj) differentiated by 
y alone and obtain 



C1D4 



1 - 



(A47) 



From here: 
1 

g'2 ~ 

^Diy^ + (Ci£>3 - C3D1) yz + CiDjy + DiGeix, z) 
Diy + {C1D3 - C3D1) z + CiDi 



(A48) 



Aft er su bstituting (|X45|) . (IX46| . ((X48)) . (|X23t and (1X24)) 
in (jA43[) we obtain 



{z - Zq) Giix, z) + {DiXo + D3Z0 + D4) — 

K 



{x - xof + [CiXo + C3Zo)^ + {z~ zaf 



+ {DiXo + D3Z + Di) Gsix, z) - ^6(2;, z) = 0. (A49) 

Following further exactly the same scheme that we pre- 
sented in solving (IA7[) we arrive at an equation determin- 
ing R{t) that is identical with (IA39[) . i.e. does not put 
any additional limitation on the R determined by (jA39p . 
Thus, (|A39p gives the final condition for the existence of 
RLPs in Case 1.1.1.1.2.1. 

Putting together all the partial results we get 



(A50) 



CM + y^- z^) 

+D3xy - Diyz -f C^D^x + {C4D3 + C5D1) y 
-CiDiZ + Di {C1C5 + C3C4) 
-EiiCiD^^CsDi), (A51) 

where U2 is given in (jA42|) . Now the equations 



dx 



Ti{x,y,z), 



dy 



Gi {x,y, z) 



(A52) 



do determine the RLPs, with J-i{x,y, z) and Giix^y^z) 
given by (IX40| - (|X42)) and (jMO t - (K5T|)-^ 

By comparing (|I54l) and (|M9| w ith (I7J5)) - (frTC|l we 
see that the latter is contained in (jA39|) as the subcase 
xo = = C3, with £)3, Di, C5 and Ei expressed in terms 
of C, Ai,A2,Bi and S2. 

Likewise, comparing (|A24I) and (|A39I) with (I9.19P - 
(I9.20[) we see that the latter is contained in (|A39I) as the 
subcase xq — Q — C3 = C^, with y^ transformed to zq 
by the coordinate transformation (y, z) = {z',y'). The 
other constants in (|A39p are expressed in terms of those 

from mim - (lon)) . 

These two special Stephani solutions will be contained 
also in several other subcases of { (|A24p . (jA39p } that will 



1 



1 



C1D3 

Di 



The correctness of l|A39|l - l|A42|l and l|A50| l - IIA5H was verified 
by the computer algebra program Ortocartan 1121 . 113(1 . 
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appear in the following text. The spherically symmetric 
Stephani solution will appear in several places as well, 
but each time it will be easy to recognise. 
For completeness, we will still consider: 

Case 1.1.1. 1.2. 1.2: ^2 = 

With Di 7^ ^ xqj and {zoj/xo,t) ,t ^ still applying, 
G2 — quickly leads to a contradiction; it is enough to 
divide (IA8|) by 2K, then differentiate the result by y and 
t and use the conclusion in the derivative by y alone. 
This case is thus empty. 

Case 1.1.1.1.2.1 is thereby exhausted, and wc go back 
to (|A28P to consider: 

Case 1.1.1.1.2.2: Di ^ 

By applying the same scheme as with Di 7^ we obtain 
formulae for J^i, Qi and R that are the subcases Di — 
of (fA39ll . (IA40|) - (IA42]) and (lA50ll - (lA5T|) . 

When Ci — 0, we obtain xo^t = 0, which contradicts 
the definition of Case 1.1.1. 

This exhausts Case 1.1.1.1, so we go back to (jA12[) to 
consider 



After substituting (jA55|) - (|A59|) in 
1 



— T5{x,z) + 
C3 2KH 

1 r. 



9j) we obtain 

xqT5{x,z) -Tq{x,z) 



{x - xo) +yo + (z - csXqY 



H — (dixo + yo + dA FAx, z) ~ 0. 

Cl 



(A60) 



Following the same sequence of steps that we described 
below (jA8l) . we arrive at: 



^ {diXQ + yo + ^4) + 1 [(1 + C32) xo^ + 



ciR ' " ' 2 

+C5X0 + ceyo = El = constant 



yo 



(A61) 



Case 1.1.1.2: {zo^t/xo,t) ,t = 
Then (|A12p becomes 



This is equivalent to the subcase C3 = of (|A39|) trans- 
formed by {y,z) = {z',y') (this transformation inter- 
changes the names of yQ and zq). Note from (|A5p that 
the interchange of y and z implies the interchange of Gy 
and Gz, and consequently the transformation {J-i, Gi) — >■ 
Th us the :F i and gi f or (|A6l|) are ob- 
tained from (|A40P - (|A42|1 and (|A50P - (jA51|l in this 
way, with the substitutions C3 — and {y,z) -> {z,y). 
To have a consistent naming of the constants, one should 
do the following replacement in (jA61l) : 



Case 1.1.1.2.1: [(l/xo^t)(l/2j^),t ] ,t ^ 

For the sake of comparison with the previously consid- 
ered cases, from this place up to (jA65l) we replace the 
capital letters denoting constants by their corresponding 
I.e. letters. Then the solution of (|A53|) is 



iDi,Ds,Di,Ci,Ci,-C5). 

(A62) 

This is the end of Case 1.1.1.2.1.1.1. We go back to 
(jASSp and consider now 

Case 1.1.1.2.1.1.2: C3 ^ 



yo,t 

Xo,t 



-(- 



,t — 



T2 



-=ci. (A54) 



Case 1.1.1.2.1.1: ci 7^ 

We then get 



Even though the limit C3 = is singular in (jA59[) - 
(|A60I) . by going through the usual procedure we end up 
with the subcase C3 = of (jA61l) . This is also a regular 
subcase of the current T\ and calculated in the way 
explained above. 

Thus we go back to (|A54p and consider now 

Case 1.1.1.2.1.2: ci = 



7^ = dxXQ + yo + 1^4, 
zA 

= (ci2/ + T^[x, z)) T2. 
Zq = C3XQ, 



(A55) 

(A56) 
(A57) 



Then ()A55P does not exist because (|A54p only deter- 
mines {yo,t/xQ.t),t = and {J^zl^2),y = 0. From the first 
of these we have 



2/0 = ^22:0, 



(A63) 



the last one from the definition of Case 1.1.1.2. 

We substitute (|A55P - (|A57P in (|Allj) . In Case 1.1.1, 
where xo.t ^ 0, the resulting equation integrates as 

c3.F1 = l + (diy- J'5(x,z))J-2. (A58) 
Case 1.1.1.2.1.1.1: C3 7^ 



and (|A56P stiU holds with ci = 0. 
Case 1.1.1.2.1.2.1: C3 7^ 

Proceeding as before by consecutive differentiations we 
end up with the following equation: 



1 



Ca 1 

+ 7^ + C2^ + C3^) XQ^ + C5XQ -f Co = 0. 



We then use (|A55p - (IA58P in differentiated by y 
and obtain 



1 _\y'^ + ^yz + d4,y + FQ{x,z) 



T2 



X — zjc-j, 



(A59) 



2KR 2K 2 

(A64) 

The function K{t) is still arbitrary at this point. Since 
in the acceleration- free limit xo and KR = fc/4 are con- 
stant, we see that k/R — constant in this limit. So, (IA64P 
can reproduce the fc = FLRW limit. 
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Running this solution through (jASp with Q2 7^ 
0, with the current values of uq and Zq, and with 
[{l/xQ^t){^/2K),t] ,t7^ as appropriate for the current 
Case 1.1.1.2.1, we obtain 



Then (lAlip becomes 



zo.t 



^2 



1 

2K 



^2 



0. 



(A70) 



1 



2KR ■ 2I + ^ + + ''"'^ + ^""^^ + ^3 = 0, 

(A65) 

where i = 1,2, 3, are new arbitrary constants. Choos- 
ing (E'l, i?2, -E3) = (c4,C5,C6) we make (|A65p identical to 
(jA6ip . Then (|A64|) does define a subcase of the Stephani 
model, different from (jA39|) . that also has RLPs. The 
corresponding J^i and Qi are found to be 



+ C2XJ/ + C3XZ + C5X + Cq 



iCg (-x2 - y2 + ^2) + + c^y^ + + C-aCe ' 

(A66) 

|C2 {x^ -y^ + z^) -xy- C^yz - C^y - CzCg 
C3 (a;2 + y2 _ ^2) _ _ c-syz - C5Z - C3C6 ' 

(A67) 

Even though the subcase C3 = requires separate 
treatment at intermediate stages of the calculation, the 
final formulae for 1/R, J-\ and Q\ turn out to be con- 
tained in (jA64[) . (jA66[) and (jA67|) as the regular subcase 
C3 = 0. While considering C3 — 0, the subcase C2 = 
also requires separate treatment, but in the end leads to 
R, J"i and given by (|A64l) ~ (|A67l) with C3 = C2 = 0. 
When C2 = 0, the Q2 in (|A8| must be zero.*^ 

When we run the general expression for (|A64p through 
(|A8P with t/2 = 0, we quickly obtain C3 = C2 = as a 
necessary consequence. 

Having thus exhausted Case 1.1.1.2.1 we go back to 
(jA53P and consider 

Case 1.1.1.2.2: [(l/a;o,t)(l ] ,t = « 



Then, using (|A53|) and the definitions of Cases 1.1.1.2 
and 1.1.1.2.2 we obtain: 



1 
2^ 



DxXq^D2, 2/0 = ^2X0, ZQ^C-iXQ. 



(A68) 

Proceeding further by the ordinary scheme we obtain 
from JH: 



(i^ixo + 2?2) + ^ (1 + C2' + C3') xo' + = C5. 

(A69) 



This is formally a subcase of (jA39[) . but it has Gx — Gy = 
Gz and so belongs in Case 2 considered further on. 

This exhausts Case 1.1.1, so we go back to (|A11[) and 
consider 

Case 1.1.2: XQ,t = 



The correctness of |)A64|| - IIA67II w as verified by the computer 
algebra program Ortocartan 



Case 1.1.2.1: j/o,t ^ 



Then we divide (jA70p by ?/o,t and differentiate the re- 
sult by t, obtaining 



yo,t 



T2 



1 

Vo.t 



1 
2K 



^3 
^2 



(A71) 



Case 1.1.2.1.1: (zo.t/yo,*) ,t 



Then we divide (|A7ip by (zo,t/2/o,t) ,t and differentiate 
the result by t. 

Case 1.1.2.1.1.1: xiW ^ 



where 



dcf 



1 

yo,t 



1 

2K 



ZQ,t 

yo,t 



(A72) 



Then {J-3/J-2),y= and (jA71|) immediately implies 
{J-1/J-2) 0. However, then (jA70p gives a contradic- 
tion with j/o.t 7^ 0- This case is thus empty, so we proceed 
to consider the complementary 

Case 1.1.2.1.1.2: xiW = 



Then we have 



1 
2K 



Di + -D22/0 + £'3^0, 



and from (|A7ip 



^ = +Ti{x,z). 

■r2 ■r2 



(A73) 



(A74) 



Proceeding from 
1 



9]) by the usual method we arrive at 



R 



-Ciyo 



D2V0 + -^32:0) 

-Cr^Z^^El. 



1 



(A75) 



This is equivalent to the subcase Ci = C3 = of (jA39p 
under the coordinate transformation {x, y) = {y' , x') that 
interchanges the names of Xq and j/o- Note, by looking at 
(|A7P - (|A8P , that the interchange of x and y implies the 
interchange of !Fi and Qi. Thus, the corresponding Ti 
and Gi are found from (IXiO]) - ([Il2|) and (|AKo1) - (|X5T|) 
as, respectively, the old Qi and J-i in the limit Ci = 
C3 = with {x,y) — > [y' ,x'), and with the additional 
renaming {C^,Di,D4) (— C5, D2, -Di) (compare (|A39p 

with (|A75jl ). 

So we go back to ()A7ip to consider 

Case 1.1.2.1.2: izo.t/yo,t) ,t ^ 
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Then 



Case 1.1.2.2: yo,t = 



zo = C'32/o- 



Case 1.1.2.1.2.1: [(l/yo,t) (l/2i^) ,t ] ,t 



(A76) 



Then (|A71|) gives 



and (|A70P imphes 



{Ml) 



(A78) 



Note that C3 7^ 0, or else (jA78p is a contradiction. 
Proceeding from (|A9p by the usual routine we obtain 



^ + § + ^ (1 + 2/0' + Csyo + = 0, (A79) 

with K{t) undetermined. This is equivalent to the sub- 
case C2 = of (IA64p . The transformation from ()A64p 
and (|A66|) - (|A67p to the current case is [x,y) = (j/', x'), 
with the accompanying renaming (a;o,?/o) — ^ (2/0, a^o)- As 
explained under (jA75p . the interchange of x and y im- 
plies the interchange of Ti and tJi, so the current JFi and 
tJi are obtained as the Qi and Fi, respectively, of (|A66p 
- (|A67P with {x, y) -» {y,x ) and C2 = 0. 

So, we go back to (jA7ip once more and consider the 
second possibility: 

Case 1.1.2.1.2.2: [(l/yo,t) (l/2j^) ,t ] ,t = 

Then 



Then = by a transformation of y, and in (IA70p 
we consider 

Case 1.1.2.2.1: z^.t ^ 



Then (IA70p may be written as 



1 



1 
2K 



0. 



(A83) 



When {(1/00,*) [l/(2/f )] ,t } ,t ^ 0, this yields {T^jT, 



2) ,y 



,y = 0, and the usual routine leads from (|A9 



to 



2KR 



2K 



+ -^ + ^^0^ + C5Z0 + Ca 



0, 



(A84) 



with K{t) undetermined. This is the subcase 03 = 6*2 = 
of (|A64p obtained by the coordinate transformation 
(x, z) = (z',x') that interchanges xq with zq. Note from 
(jASP that the interchange of x and z implies the inter- 
change of Gx and Gz, and consequently the transforma- 
tion {J-i,Gi) i^/ J-iiGi/ J-i)- Thus, the J^i and Qi for 
(jA84p are obtained from (|A66p - (|A67p in this way, with 
the substitutions C2 = C3 = and {x, z) — > (z, x). 



Thus we assume 



^7(27?)'* -*=Oin M. Then 



2K = ^1+^3.0, 

^ = + T4:{x,z). 

J-2 J-2 



(A85) 
(A86) 



— = Di + D22/0, 



and from (|A70I) . since yo,t 7^ 0, 



= £'2^ - y + ^"4(2;, z). 
"''2 -^2 



(A80) 



(A81) 



Case 1.1.2.1.2.2.1: C3 7^ 

Then, by the usual routine, (IA8ip used in (|A9p leads 

to 

^ + 1)22/0) + ^ (1 + C^) yo^ + Ciyo = C5, (A82) 

which is formally the subcase C5 = -D3 = 0, zo = Cs^o of 
(IA75p . However, in this case we have Gx = Gy = Gz = 
in (|11.5p - (lll.6p . Thus, it is in fact excluded from 
the present consideration and will appear later, when we 
consider Case 2. 

Case 1.1.2.1.2.2.2: C3 = 



Then the subcase C3 = of (|A82I) is obtained, so it 
belongs in Case 2, too. 

This exhausts Case 1.1.2.1, so we go back to (jA70p to 
consider 



This leads from (lA9l) to 



R 



{Di + A3^o) 



-zo2 + C5Zo = £;i. (A87) 



This is the subcase yo = of (|A75p . As with (|A82p . 
in this case we have Gx — Gy = Gz = 0, and so it is 
excluded from the present consideration - it will appear 
in Case 2. 

So we go back to (jA7ip and consider 

Case 1.1.2.2.2: zq * = 

With xo^t = yoj — zqj — now being the case, we 
are in the spherically symmetric subcase of the Stephani 
solution. Then (jA9p may be written as 



X 



>2 



+ - ix:' + + z^) + 

2KR 2 ^ ^ ' 2K F2 



The t-derivative of this is 



/ 1 



\2KR 



(A88) 



(A89) 



If {2KR),t = 0, then this is simply the FLRW family of 
models that we need not investigate. Consequently, we 
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take {2KR),t 7^ 0. Then both the other factors must be 
nonzero, and 



1 



2KR 



This is the subcase of the sphericahy symmetric 
solution that we identified in Sec. |Vl and the 
above is consistent with (|5.10p . AU of its nuU 
are RLPs. The FLRW hmit of this model is Ci 
it includes the whole FLRW family. 

Equation (|A90|) is the final solution of (|X7t . 
vent (IA8|) from imposing any limitations on it, 
ficient to choose G2 — G3 — 0, Qi — y/z. 

This exhausts Case 1.1, so we go back to (|A7p 

Case 1.2: ^2 = 

Then (IA7I) can be rewritten as 



(A90) 

Stephani 
equation 
geodesies 
= 0, and 

. To pre- 
it is suf- 

and take 



X ~ Xo ^ {z ~ Zo) -Fi + 7777-^3- 
ZK 

We differentiate this by y and t and obtain 
- zo.tJ^i.y + ( ^ ) '* -^a.iy = 0. 



Case 1.2.1: zo,t ^ 
Then from (|X9T|) 



(A91) 



(A92) 



(A93) 
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: 7^ 0, we have Ti^y = T^^y = 0, and 



When 

then we differentiate (|A9ip by x to get 

1 



1 = (z - Zo)Ti^a: ■ 



2K 



3, a; ■ 



(A94) 



Differentiating this by t we get Ti^x = ^3,x = in con- 
sequence of the assumptions about the functions of t, 
but this is a contradiction with (IA94p . This means that 



zo.t (2_R-) 



0, i.e. 



— D'^Zn 

2K ^ " 



Da 



(A95) 



Continuing from (IA91[) by the usual procedure we get the 
following final solution; 



2^0 — C2Zq, 

D3X + C2Di 



Ti = 
^3 - 



D4 + D3Z ' 

X — C2Z 



D4 + D3Z' 

and we next have to verify (jA8[) . We first try 
Case 1.2.1.1: ^2 7^ 



(A96) 
(A97) 

(A98) 



Then 

y-yo 



can be written as 
1 



G2 =^''''^¥2 



2KR 



1 r 



{x ~ C2Zaf + (y - yaf + {z~ z^f 



(A99) 
2K ^2' 



Using (IA95|) and (|A96p we take the second derivative of 
this by y and t. Since we are in Case 1.2.1, where zq^j 7^ 0, 
we get 



yo.t 

ZQ,t 



1 



-1 



(AlOO) 



Case 1.2.1.1.1: {yo^^/ zoj) ,t ^ 



Then (|AlOOp solves as 

y + g4{x,z), 
Gi = D3Q3 



1 

G'2 



Q5{x,z) 



y + Qi{x,z)' 



(AlOl) 

(A102) 



Proceeding from (jA99|l by the usual routine we find ^^3 
and the arbitrary functions of (x, z), and we end up with 



{Di + D^z^) + i [2/0' + (1 + C2') zo^ 
+C42/0 + Cszo = El. 



(A103) 



This is equivalent to the subcase I?i = Ci = of (jA39|) 
under the coordinate transformation (x, y) — (y', x') that 
interchanges xq with j/q. The remark under (jA75p about 
the transformation of J-i and Qi applies also here. 
We go back to (jAlOOp and consider 

Case 1.2.1.1.2: (j/o.t/zp^t) ,t ^ 

Then 

2/0 = C3Z0, (A104) 

g^ix.z). (A105) 

In the usual way this leads from (IA99P to 

4 [Di + D3Z0) + i (1 + C2' + zo^ - C5Z0 = Si, 
H Z 

(A106) 

which is formally the subcase Ci — C4 — Di ^ 0, xq — 
C2Z0 of (|A39p . but belongs to the class with = Gy — 
Gz = considered under Case 2. 

This exhausts Case 1.2.1.1., so we consider 



Qi ^ Ds^ + Ca (—-y 
Q2 ' G2 \G2 



Case 1.2.1.2: ^2 







Then (K8\i becomes, using (|A95P 

y-y„^{z- Zo) Gi + {D4 + D3ZQ) G3- (A107) 

By the normal routine this leads to 

yo = C3Z0, (A108) 
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Qi = 



Day + C^Dj 
Di + D3Z 

y ~ C3Z 

Di + D3z' 



(A109) 
(Alio) 



which does not impose any limitation on (1X95]) - (|X98)) 
and leaves R undetermined. Thus, this subcase has a 
nontrivial FLRW hmit. We have 



V = 



1 
R 



+ —[x^+y^ + z'^-2 {C2X + Csy + z) zq 



4i? 



2R 



D4 + D3Z0 



(Alll) 
(A112) 



The flat FLRW model is contained here in the limit 
(£'4, D2) — 00. In the general case, an orthogonal trans- 
formation of (x, y, z) may be used to achieve C2 = C3 = 

{zq is then transformed to 29 



+ Cg^zo), and 



then the model is seen to be axially symmetric, with the 
orbits of symmetry in the new (x, y) plane. Assuming 
C2 = C3 = and using (|M09l) and (PlTT|) in 

(|XT|) we get 



d^ 



D3X 



D4 



D.z' 



dy 
dz 



Day 



D4 + D3Z 



(A113) 



as the equations defining the RLPs in this case. As is 
easy to see, they obey ydx/dz — xdy/dz — 0, and so are, 
in the {x, y) surface, straight lines passing through the 
symmetry axis x = y = 0. 

Thereby, Case 1.2.1 is exhausted, so we go back to 
(|A92|) to consider 

Case 1.2.2: zo,t = 



This is equivalent to zq = 0. Then (IA92[) leads to two 
further subcases: 

Case 1.2.2.1: (1/2X) ,t 7^ 



By the usual method we obtain from (jA92|) and (jA91|) : 

- - E4 



-^3 = £^3, 



E3_ 

2K 



E4-X0, (A114) 



with R undetermined. This is the final solution of (jA7p . 
For checking (jA8[) we have to consider separately 

Case 1.2.2.1.1: £3^0 



Then (|All4p determines K, and from the definition of 
Case 1.2.2.1 it follows that a;o,t ^ 0. But to continue, we 
have to separately consider Q2 being zero or not. 

Case 1.2.2.1.1.1: ^2 7^ 

Then (|A8|) is written as 



y-yo 



Q2 
1 

'2 



2KR 



{x - xo) + [y - yo) + z + — (£4 - a;o) 



(A115) 

^3 



Taking the derivative of this by y and t, and knowing 
that XQ^t 7^ 0, we get 



yo,t 
xo,t 



^ ,„-l 



£3 \Q2 



(A116) 



Case 1.2.2.1.1.1.1: {yo,t/xo,t) 



Then 
1 



y + G4{x,z), 



G3 



G5{x,z). 



(A117) 



G2 G2 
Proceeding from (jAllSp in the usual way we obtain 

-J— {E4 ~ Xo) + 1; (xo^ + yo^) + C4X0 - Csyo = Ei. 
j^^K z 

(A118) 

This is the subcase of (jA39|) that results when we take 
Ci = C3 = -D3 = in (jA39p and interchange yo with 
zo- How the new J-i and Qi are calculated after such a 
transformation is explained under (IA6ip . To have a con- 
sistent naming of the constants, one must take in (|A39p 
{Di,D4) = {-l/ E3 E^/ Es). 
We go back to (|A116p and consider 

Case 1.2.2.1.1.1.2: (yo,t/xo,t) ,t = 



Then from the definition above 
yo = C2X0, 
and by integrating (|A116P 

^3 



— C2E3 



1 
G'2 



y 



EiG4{x, z). 



(A119) 



(A120) 



By continuing from (jAllsp in the usual way we get finally 
{E4 -xo) + \{l + C2^) xo^ + C4X0 = El. (A121) 



This is formally the subcase yo — C2X0 of (jA118p . again 
with some renaming of the constants. However, it has 
Gx = Gy = = 0, so in fact it belongs in Case 2 
considered further on. 

This completes Case 1.2.2.1.1.1, so we go back to 
(|A114p to consider 

Case 1.2.2.1.1.2: ^2 = 



Then (|A8P is written as 
1 



2/0 = 2:^1 + 7;- {E4 - a;o) Ga- 
E3 



(A122) 



Case 1.2.2.1.1.2.1: xo,t ^ 

Then the final solution of (|Al22p is 



yo = C2X0, 



Gi = -{y-E4C2) 



G3 — C2E3, 
(A123) 
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where C2 is allowed to be zero. This imposes no ad- 
ditional conditions on (|A114p . and R remains undeter- 
mined. The function V is in this case 



1 



2 I 2 I 2 

X + y + z 



2 {Ei - xo) 
2x^ [x + C^y) + (1 + 



This is equivalent to (jAlll[) - (|A112p , which is seen when 
we transform C2 to by a rotation in the {x, y) plane. 
Subsequently, the new (x, z) have to be transformed by 
(a;, z) = (z',a;'). How the new Fi and Qi are calculated 
after such a transformation is explained under (jA84[) . 
We go back to (jA122p to consider 

Case 1.2.2.1.1.2.2: t = 



Then (|A122p implies ?/o,t = 0. With xq,* = yo.t = 
zo,t = we are then in the spherically symmetric subcase 
of the Stephani solution, but not in its full generality. 
Since (IA114I) still applies, with xq = constant we have 
K = k/{AR) ~ constant, and this corresponds to the 
subcase A2 = of (|5.10|) . Since this has already been 
investigated, we go back to (jA114[) to consider 

Case 1.2.2.1.2: E-i = Q 



Then (|A114p leaves K still arbitrary, but implies xq 
Ej^. By a transformation of x this can be changed to 



xq — Ei — 0. 



(A125) 



With this, we go on to consider (|A8|) . Recall that we are 
still in Case 1.2.2, where also zq = 0. 

Case 1.2.2.1.2.1: ^2 7^ 



Then 



is written as 



y-yo 



+ - 



Q2 
1 



2KR 



(A126) 



+ {y- yo) + 



j_g3 

2K 02 ' 



Taking the derivative of this by y and t we get 
1 



-yo,t 



(A127) 



We are still in Case 1.2.2.1, where {1/2K) ,t ^ 0, so we 
divide (|A127p by {1/2K) ,t and consider first 



Case 1.2.2.1.2.1.1: [yo.t/ {1/2K) ,t] ,t ^ 



2) ,y 



0. 



Then (|A127|) implies (l/^2),,j-l 

This leads to a simple coordinate transform of (jA84[) [zq 
replaced by yo)j which is a subcase of (|A64| . 

We go back to (|A127|) and consider 

Case 1.2.2.1.2.1.2: [yp^t/ {1/2K) ,t ] ,i = 

We have then j/q — E)i/{2K) + D2- We do not have to 
consider the case Di = because this would mean y^ — Q 
in addition to xq = zq = 0, and we would be back in the 



spherically symmetric Stephani solution, considered in 
Sec. |Vl So the new situation arises only when Di ^ 0, 
and then we can rewrite the last formula as 



1 

— = Di+D2yo, 



(A124) obtaining from (KT27\ 



7^ ^ y - E>2^ +Gi{x,z). 

b'2 t'2 



(A128) 



(A129) 



This leads to the subcase yo = of (|A118p . with xq 
subsequently changed to yQ. This again has = Gy = 
Gz — and belongs in Case 2. 

Case 1.2.2.1.2.1 is now exhausted, so we go back to 
(jA125[) and consider 

Case 1.2.2.1.2.2: G2 ^ 



Then (IMl) becomes 



This leads to 
Gi = , 



(A130) 



^3=C^5, ^=^4-yo. (A131) 



When C5 7^ this is equivalent to (IA123I) - (|A124|) . with 
the constants in (jAlSlI) related to those in (jA123|) by 
(C4, C5) = C2{Ei, E3). When C5 = 0, we are back in the 
spherically symmetric subclass. 

This completes Case 1.2.2.1. We have to go back as 
far as (jA92[) and consider 

Case 1.2.2.2: (l/2X),t= 



With zo,t = now being considered (Case 1.2.2), (IA92[) 
is fulfilled identically, and (IA9ip immediately implies 
xo.t ~ 0. Taking xq = zq = and 1/2K ~ B ~ con- 
stant in (jA9ip we get 



X — zT\ + BTz , 



(A132) 



with no limitation on i?, and this is the final solution of 
[TJ. Running this through (IA8|) with CJ2 7^ we obtain 



^ + \y^^ + ^42/0 = 



(A133) 



This case again has Gx — Gy — Gz — 0, and so belongs 
in Case 2. 

Running (|A132I) through with G2 ^ ^ we imme- 

diately get yo,t = and no limitation on R. This is a 
subcase of the spherically symmetric Stephani solution, 
the same one that we obtained in Case 1.2.2.1.1.2.2. 

This completes Case 1.2, and the whole Case 1. So we 
go back to (Ill.Sp - (|11.6p and consider 

Case 2: = 



Then we immediately have Gj, = Gj^ = in (|11.5I) - 
(dUl). This means that (|TT3)) and (pT^ are fulfilled 
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for any dx/dz and dy/dz, i.e. that all null geodesies are 
RLPs in this case. 

Multiplying the equations Gx ~ Gy ^ Gz ^ hy V,^ 
we turn them into polynomials in {x,y,z). Taking the 
coefficients of x'^ in each polynomial we get^ 

-2KK,u XQ^t + 2Kxo,ttK,t +4Kjxo^t = 0, (A134) 
2KK,tt yoj - 2KyoMK,t -4.K/yo^t = 0, (A135) 
2KK,u zo,t - 2Kzo,ttK,t ~4.K/zoj = 0. (A136) 

These are easily integrated, but a few cases have to be 
considered separately. One solution is xo,t = yo.t = ^o,t = 
0, but this is the spherically symmetric Stephani model 
that we investigated in Sec. |Vl So we assume that at 
least one of the functions {xq, yoi zq) is non-constant. By 
a coordinate transformation we may choose this to be xq. 
When xo,t 7^ 0, (IA134|) is integrated with the result 



+2K^xo.t {xo/ + yo/ + zo/) = 0, 



(A140) 



and the analogous equations with xq,* replaced by t/o,t 
and zq t ■ 



Using (Ul37l) and ^MM, eqs. (fUSOl) and (|Il40l) 
become identical and are integrated with the result: 



^ {DiXo + D2) + (1 + + Xo' + ^4X0 - C5. 

(A141) 

By an orthogonal coordinate transformation in the 
(x,?/, 2) space one can obtain C2 ^ C3 ~ 0. Then the 
model is seen to be axially symmetric, but is more general 
than the plane- or hyperbolically symmetric subcases. 
Case 2.2: Di ^ 



= D1X0 + D2, 

and then we consider 
Case 2.1: Di ^ 
In this case, (|A135P - (|A136p imply 

yo = C2X0, zo = C3X0, 



(A137) 



(A138) 



with zero values of C2 and C3 allowed. 

The terms free of {x, y, z) in the equations Gx — Gy 
Gz = are all the same: 



Then_^= 1/(202) = constant, and (jA134|l - (|A136|) 
and (|A139P are fulfilled identically. Instead of (|A140j) 
and its associated equations we then have their subcases 
K^t = which imply either xqj — ?/o,t = zo^t ~ or 
(|A138P plus the integral of (IA140p . which is 



4 + ^ (1 + ^2' + xo^ -I- E^XQ = E2. (A142) 
H IU2 



This is the subcase Di = of (|A141|) . 



-2KK,t [xo/ + y^/ + 20,(2) = 0. (A139) 



Likewise, the terms linear in {x^y, z) are the same in all 
3 equations, namely: 
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